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As is well known, approximate integral equations for liquids, such as the hypernetted chain (HNC) 
and Percus- Yevick (PY) theories, are in general thermodynamically inconsistent in the sense that 
the macroscopic properties obtained from the spatial correlation functions depend on the route 
followed. In particular, the values of the fourth virial coefficient B4, predicted by the HNC and PY 
approximations via the virial route differ from those obtained via the compressibility route. Despite 
this, it is shown in this paper that the value of B4 obtained from the virial route in the HNC 
theory is exactly three halves the value obtained from the compressibility route in the PY theory, 
irrespective of the interaction potential (whether isotropic or not), the number of components, and 
the dimensionality of the system. This simple relationship is confirmed in one-component systems by 
analytical results for the one-dimensional penetrable-square-well model and the three-dimensional 
penetrable-sphere model, as well as by numerical results for the one-dimensional Lennard-Jones 
model, the one-dimensional Gaussian core model, and the three-dimensional square-well model. 



I. INTRODUCTION AND STATEMENT OF 
THE PROBLEM 

According to equilibrium statistical mechanics^— all 
the relevant structural and thermodynamic properties of 
a one-component fluid made of particles interacting via 
a pairwise potential </>(r) = <f>(— r) are contained in the 
pair correlation function g(r; p, j3) = g{—r; p, /3), where p 
is the number density and j3 = l/fc^Tis the inverse tem- 
perature. In particular, the thermodynamic quantities of 
the fluid can be obtained from <?(r; p, /3) through different 
routes. The most common ones are the virial route, 

Z<j>,P) = &j. = l + £ J dry(r;p,/3)r.V/(r;/3), (1.1) 

the compressibility route, 

= l + pjdr {[/(r; (3) + 1] y(v; p, (3) - lft.2) 

and the energy route 

u{p,P) = ^ - f / dryfr^A/fr;/?). (1.3) 

In Eqs. (|l.ip - (|1.3p . p is the pressure, Z is the compress- 
ibility factor, d is the dimensionality of the system, \ is 
the isothermal susceptibility, u is the internal energy per 
particle, /(r) = e _/3l ^ r ) — 1 is the Mayer function, and 
y(r) = exp[/3(f)(r)]g(r) is the so-called cavity function. 
Note that in Eq. ([TT2]) the integrand [/(r) + l]y(r) - 1 = 
g(r) — 1 = h(r) is the total correlation function. It is re- 
lated to the direct correlation function c(r) through the 
Ornstein-Zcrnike (OZ) relation^ 

h(r) = c(r) + p [ dr' h(r - r')c(r'). (1.4) 



In terms of c(r), the compressibility equation of state 
(|1.2[) can be rewritten as 

X- 1 = 1 - pc, (1.5) 

where 

c = J drc(r) (1.6) 

is the spatial integral of c(r) or, equivalently, its Fourier 
transform at zero wavenumber. 

Thermodynamic relations establish that 

X- 1 (p,P) = -0-[pZ(p,P)}, (1.7) 



p— u (p,l3) = —Z(p,l3). (1.8) 

Of course, the exact correlation function g cxllct (r; p, f3) 
yields functions Z(p,(3), x(P>P)> an d u(p,j3) satisfying 
Eqs. (jl.7|) and (|1.8[) . On the other hand, except in one- 
dimensional problems restricted to nearest-neighbor in- 
teractions, the function g exact (r; p, j3) is not known and so 
one has to resort to approximations. The price to be paid 
is that, in general, an approximate function g app (r; p, j3) 
violates Eqs. (|1.7[> and (|1.8p and so it is thermodynam- 
ically inconsistent. There exist, however, a small num- 
ber of approximations in which the energy route is con- 
sistent with the virial route, i.e., Eq. (|1.8[) is verified. 
This class of approximations includes the hypernetted 
chain 4 (HNC) and the linearized Debye-Hiickel^ theo- 
ries for arbitrary potentials, the mean-spherical approxi- 
mation for soft-potentials r- and the hard-sphere limit of 
the square-shoulder potential for any approximation. 7,8 
On the other hand, even in those cases the virial- 
compressibility consistency condition, Eq. (|1.7p . is not 
satisfied. In fact, some liquid state theories include one 
or more adjustable state-dependent parameters which 
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are tuned to achieve thermodynamic consistency be- 
tween two or more routes. This is the case, for in- 
stance, of the modified HNC closure,— >^ the Rogers- 
Young closure the Zerah-Hansen closure^ the self- 
consistent Ornstcin-Zcrnikc approximation^ the hier- 
archical reference theory,— Lee's theory based on the 
zero-separation theorems^ the generalized MSA^ or 
the rational- function approximation. 17 

The thermodynamic inconsistency problem manifests 
itself at the level of the virial coefficients B n (j3). They 
are defined by the series expansion 



1 d 



]() pB n M n>2. 



(1.18) 



Z{p,p) = l + Y^B n+l {P)p n . 



(1.9) 



Likewise, the series expansion of the isothermal suscep- 
tibility, the internal energy, and the cavity function in 
powers of density can be written as 



n=l 



P n=l 



y(r;p,/3) = l + ^y n (r;/3)p". 

n=l 



(1.10) 



(1.11) 



(1.12) 



Here the subscripts "c" and "e" denote virial coefficients 
obtained via the compressibility and energy routes, re- 
spectively. Obviously, the exact functions y® xact give con- 
sistent virial coefficients, i.e., B™ ct = B^ ct = B^ ct . 
As said above, the HNC theory provides consistent ther- 
modynamic properties through the energy and virial 
routes only, i.e., B™C = gHNC ^ b hnc Except for 

a few cases,— ~— an approximate theory generally predicts 
different sets of thermodynamic quantities from the three 
routes. This is the case of the well-known Percus-Yevick 
(PY) theory^ i.e., B^ ? i?£ + B™. 

Most of the liquid state theories, including HNC and 
PY, give the exact cavity function to first order in 
density and hence provide the exact second and third 
virial coefficients. Therefore, the fourth virial coeffi- 
cient is usually the earliest one revealing the approx- 
imate nature of the theory. In principle, the two 
HNC coefficients (Bf® c = Bf® c and B^ G ) and the 
three PY coefficients (B^, B^, and B^) are un- 
related. On the other hand, recent analytical eval- 
uations of the fourth virial coefficient for the three- 
dimensional penetrable-sphere potential^ and for the 
one-dimensional penetrable-square-well potential^ show 
a remarkably simple relationship between B™ C an d 



£?4 (f, namely 



Inserting Eq. (|1.12p into Eq. (11.11) one can express B n in 
terms of an integral involving y n -2, 



(1.19) 



£f n ,vC8) = — / dr y n - 2 {r; f3)r ■ V/(r; (3), n>2, where, according to Eq. (11.170, 

Zd J 



(1.13) 

where we use the subscript "v" to emphasize that Eq. 
(|1.13[) gives the nth virial coefficient through the virial 
route. Similarly, insertion of Eq. (|1.12[) into Eqs. (II. 2[) 
and (|1.3|) yields 



B, 



4 (2X2X3 ~ xl 



Xi) 



(1.20) 



X aC8) = / dr/(r;/3), 



(1.14) 



Xn(P)= / dr [f(r;P) + l]y n - 2 {r;P), n > 3, (1.15) 



u n (0) = -\f dr t/„_ 2 (r; 0)^/ (r; /3), n > 2. (1.16) 



Use of the thermodynamic relations (jl.7l) and (11.8[) al- 
lows one to express the virial coefficients B n from the 
coefficients \n an d u n , respectively, as 

^ n— 1 

B n ,M = J2 mB m , c {P) Xn +i- m {P), n > 2, 

n * — ' 

m— 1 

(1.17) 



The two potentials considered in Refs. I18U19I have two 
common features: they are (i) bounded and (ii) stepwise 
constant. Thus, it might be the case that the relation 
()1.19j) is closely connected to one or both features, not 
being generally valid. The aim of this paper is to show 
that this is not the case by proving Eq. (|1.19[) for any 
dimensionality d and any potential <f>(r). Furthermore, 
we will show that Eq. f|l . 19[) can be extended to fluid 
mixtures with any number of components. 

This paper is organized as follows. The diagrammatic 
representation of the cavity function to second order in 
density and of the fourth virial coefficient is presented in 
Sec. HIl This is followed by the mathematical proof of 
Eq. (|1.19l) , where use is made of identities derived in the 
Appendix. Section IIIII provides a few examples where 
Eq. (I1.19p is numerically verified. Next, it is shown in 
Sec. II VI that the results remain valid in the more general 
case of a multicomponent fluid. Finally, the results are 
discussed in Sec. IVl 



II. PROOF OF THE RELATIONSHIP 



B. 



HNC 



3d 



PY 



Substitution of Eq. (p3]l into Eq. (fTTT3|) for n 

yields 



Statistical-mechanical methods allow one to express 
the functions y n (r; P) as sums of n-particle multiple inte- 
grals of products of Mayer functions. These integrals are 
conveniently represented by diagrams.— ~— In particular, 



yi(r;/3) 



y exact (r;/3) = _ 2 




(2.1) 




(2.2) 



The open circles represent two fixed root points separated 
by a vector ri2 = r, the filled circles represent field points 
to be integrated out, and each bond represents a Mayer 
function. For instance, 




= / d3/i 3 / : 



2.S. 




d3 / d4/ 13 / 34 /24/l4, 



(2.3) 



(2.4) 



where di = dr^ and fij = fji = /(r^ — rj; f3). The factors 
2 and 4 affecting the first and second diagrams of Eq. 
(|2.2[) . respectively, reflect the degeneracies of those dia- 
grams with respect to the exchange of field points or of 
root points. As will be seen in Sec. lIVl these degeneracies 
are broken down in the multicomponent case. 

Both the HNC and the PY approximations are consis- 
tent with Eq. (I2TT1) but not with Eq. (12T2")) . In the case 
of HNC, the last diagram on the right-hand side of Eq. 
(|2.2[) is neglected, while PY neglects the last two dia- 
grams. More in general, imagine an approximation that 
includes those last two diagrams but with weights Ai and 
A2, respectively: 



y { 2 XlM \r;(3) 




The class of functions y^ 1 '^ (r; /?) includes the exact, 
HNC, and PY functions as special cases: 



2/2 Xact (r;/?) 

y 2 HNC (r;/3) 
2/P>;/3) 



yr\r;f3). 



(2.6) 



Bt M \P) = Td 2 



+Ai 




(2.7) 



where here one of the root points has become a field point 
and a dashed bond represents a term of the form r-V/(r). 
For instance, 



d2 / d3 / d4/ 13 /34/24/l4ri 2 -Vl2/l 




(2.8) 

Note that the one-root diagrams do not depend on the 
location of the root, so they can be expressed as zero- 
root diagrams divided by the d-dimensional volume V. 
For example, 



(2.9) 



In the case of the isothermal susceptibility coefficients, 
Eqs. Q| , (Ti~T5l) . (EH), and give 






X2 = 



(2.10) 



X3 = 




In Eq. (|2.12[) we have taken into account that 





(2.13) 



Equations (fi~20| and (f2~T0|) - (|2~T2|) readily give 



B ( t M) {P) 



2 + Ai 



4 + Ai + A 2 



where use has been made of the properties 



X2X3 




xl 



(2.14) 



(2.15) 
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Equation (|2.14[) can also be obtained from Eq. (jl.5l) . This 
is the path followed in Sec. IIVI for the multicomponent 
case. 

Now, in order to connect Eqs. (I2.7[) and (|2.14[) . we 
need to make use of the following properties 2 - (see the 
Appendix for a derivation): 
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(2.16) 



(2.17) 



(2.18) 



Therefore, Eq. (|2.7I) can be rewritten as 



B 



(Ai,A 2 ) 
4,v 
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Ai 
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4rf 



(2.19) 



Comparison between Eqs. (|2.14|l and (|2.19p shows that 

3 



2 + Ai 



-B 



4,c 



03), 



provided that 



A 2 = 



3Ai 



2 + A : 



(2.20) 



(2.21) 



This is the main result of this paper. If Ai = 1 then 
A2 = 1 and we recover the exact consistency condition: 
B exact = S exact Qn the oth er ha nd, the choice X x = 
yields the sought result, Eq. (|1 . . 

It is quite apparent that the Eq. (|2.20|) is, from a math- 
ematical point of view, more general than (|1.19[) . Let 
us define the (Ai, Ai)-subclass of approximations as the 
one compatible with Eq. (|2.5I) with < Ai = A2 < 1. 
Analogously, we define the (1, A2)-subclass of approxi- 
mations as the one with A x = 1 and < A 2 < 1. 
Obviously, any approximation retaining the exact 1/2 (f) 
(i.e., Ai = A2 = 1) belongs to both subclasses. On 
the other hand, the PY and HNC theories are mem- 
bers of the (Ai, Ai)-subclass and (1, A2)-subclass, respec- 
tively. Equations (|2.20p and (|2.21[) then state that, for 
any (Ai, Ai)-subclass approximation, there exists a spe- 
cific (1, A2)-subclass approximation, such that the com- 
pressibility and virial values, respectively, of -B4 are pro- 
portional each other. The connection between both sub- 
classes is schematically illustrated in Fig.[TJ Interestingly, 
the largest deviation of the proportionality factor from 1 
occurs at Ai = 0, i.e., in the case of the PY and HNC 
pair. 




FIG. 1: The diagonal (labeled c) and vertical (labeled v) 
lines represent the (Ai, Ai)-subclass and the (1, A2)-subclass 
of approximations, respectively. Both lines intersect in the 
exact case Ai = A2 = 1, while the other ends correspond to the 
PY and HNC approximations, respectively. The virial-route 
value of the fourth virial coefficient in any (1, A2)-subclass 
approximation is related to the compressibility-route value in 
a given (Ai, Ai)-subclass approximation. This is represented 
by the dots connected by dashed lines. 



III. NUMERICAL EXAMPLES 



As said in Sec. HJ the simple relationship (|1.19l) de- 
rived in the preceding section is confirmed by analyt- 
ical results for the three-dimensional penetrable-sphere 
model 1 - and for the one-dimensional penetrable-square- 
well potential. 19 The analytical derivation of the exact, 
HNC, and PY expressions of the fourth virial coefficient 
was possible in Refs. [l8lfl9l thanks to the simplicity of 
the interaction models (stepwise constant). In general, 
however, the computation of B4 is a complicated task 
that requires numerical methods, especially in the case 
of continuous potentials. 

As a numerical test of Eq. (|1.19|) , we have numerically 
evaluated Bf xact , B™ c , B™ c , SfJ, and for two 
one-dimensional interaction models. The first one is the 
conventional Lennard- Jones (LJ) potential: 



(r) = 4e 



(7)' 



(3.1) 



This is an unbounded potential, repulsive for r < 2 1 / 6 <r 
and attractive for r > 2 1 / 6 <r. As a second example, we 
have considered the bounded, purely repulsive Gaussian 
core model^ 



(r) = eexp [-(r/cr) 2 ] 



(3.2) 



The results are displayed in Figs. [H and [31 respectively. 
In the case of the LJ potential, Fig. [2] shows that the best 
general agreement with the exact results are provided by 
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FIG. 2: Temperature dependence of the fourth virial co- 
efficient for the one-dimensional LJ potential. The thin line 
with triangles represent the exact results, while B™ c , B™ c , 
BfJ, and BfJ are represented by the thick solid, dotted, 
dashed, and dash-dotted lines, respectively. The circles rep- 
resent the values of f-B^, which fall on the B™ c curve. 



j3 




kj/e 

FIG. 3: Same as in Fig. [2] but for the one-dimensional 
Gaussian potential. Note that the B™ c and B^ curves 
are hardly distinguishable in this case. 



B™ c and the worst by B^ . As the temperature in- 
creases, i?4 Y eventually becomes quite accurate. What 
is more important from the point of view of this work is 
that the numerical values of §-64 J are practically indis- 
tinguishable from those of B™ c , in agreement with Eq. 
(PI). 

Regarding the Gaussian core potential, Fig. [3] shows 
that i?4(f provides the best description for fcsT/e < 0.9, 
while B™ c does for fc^T/e > 0.9. It is also notewor- 
thy that B™ c ~ Bf y for the whole temperature range 
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FIG. 4: Same as in Fig. [2] but for the three-dimensional SW 
potential with A = 1.5. The data have been extracted from 
Fig. 10 of Ref. HI 



shown. Again, the numerical results confirm Eq. (|1.19[) . 

As a final test, we have turned our attention to the pio- 
neering numerical evaluation by Barker and Henderson 21 
of the fourth virial coefficient (they also considered the 
fifth one) for the three-dimensional square- well (SW) po- 
tential 



oo, r < a, 
(r) = < -e, a < r < Act, 
0, r > Act. 



(3-3) 



The results for a potential range A = 1.5 are displayed 
in Fig. 2J In this case, B™ an d give the most ac- 
curate results for 1 < fc^T/e < 1.5 and ksT/t > 1.5, re- 
spectively. Moreover, both approximations are observed 
to be consistent with Eq. (| 1 . 19[) within unavoidable nu- 
merical uncertainties. 



IV. EXTENSION TO FLUID MIXTURES 

So far, in Secs.lB lIlTl we have assumed a one-component 
fluid. However, the main result of this paper, Eq. (|1.19p . 
can be extended to the case of a mixture, as shown below. 

Let us consider a multicomponent fluid made of an 
arbitrary number of species with mole fractions Xi, at 
density p and inverse temperature /3. The interaction 
potential between a particle of species i and a parti- 
cle of species j is <j)ij{ r ) = <^-j(— r). The pair cor- 
relation and cavity functions are <7ij(r; p~, /3; {xk}) and 
yi](?;p;P;{x k }) = exp[f3(f> l: j(r)}g i:j (r;p;P; {x k }), respec- 
tively. The multicomponent version of Eq. (jl.ip is simply 



1 



-Y 



XiXj J dr yij(r)r- V/#(r). 



(4.1) 
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The compressibility equation of state for the multicom- 
ponent case is more conveniently expressed as a general- 
ization of Eq. (|1.5j) rather than of Eq. (|1.2[) . It is given 
by ? 

x^ 1 = i - /'X ( 4 - 2 ) 

where Cy = Cji are spatial integrals of the direct correla- 
tion functions Cy (r). The latter are defined in terms of 
the total correlation functions hij (r) = (r) — 1 via the 
OZ relation 

hj j (r) = Cjj (r) + p Xk / dr' h ik (r - r')c kj (r'). (4.3) 
k •* 

This equation yields 

= Cij + p 2J XkhikCkj, (4.4) 
fe 

where again the tilde denotes spatial integration. 

Let us now consider the virial expansion. The virial 
coefficients B n are still defined by Eq. (jl.9[) . but now they 
can be expressed in terms of composition-independent 
coefficients. For instance, 

B 2 = ^2x i XjB ij , (4.5) 

id 

B 3 = ^ XiXjXkBijk, (4.6) 

hi 

Bi = 2J XiXjXkXiB tJ k£. (4.7) 

i,j,k,t 

As for the correlation functions, the generalization of Eq. 
(fL~T2|) reads 

n=l 

A similar expansion can be carried out for ft-y(r) and 
Cy(r). In particular, 

hij-fi=Jij, (4.9) 



= J/ijjn + " > 1, (4.10) 

Cij-fl = %j, (4.11) 

Cfj;l = Vij:l + ~ J! X kf*kfkj, (4.12) 

fe 



Cy;2 = Vij;2+fVij;2-^Xkfik(ykr,l+fVkj-,l) 

k 

~^2 x kfjk (yki:l +fUki;l) + ^2%kXtfikfktftj- 

k k.e 

(4.13) 

In these equations, fUij. n denotes the spatial integral of 
fij{Y) yij , n {v) 

The coefficients yij- n {v) in Eq. (|4.8[) are polynomials of 
degree n in the mole fractions. In particular, 

yij-A r ) =^2%kyij\k(*), (4-i4) 

VijM T ) = y^, x kxeyjj\ke(r)- (4.15) 

k,t 

Analogously, 

Cy;2(r) = ^XfcX£Cy| W (r). (4.16) 

k,e 

From Eqs. (|1.7|) and (|4.2[) it follows that the fourth virial 
coefficient provided by the compressibility route is 

B 4 , c = -~ y^XjXjCjj^, (4-17) 

so that 

Bijke,c = —-^{cij-ki + Cik-ji + cn-jk + Cjk-u 

+Cji-ik +C k £-ij) . (4-18) 

The function yij\ k (r) is represented by the same di- 
agram as in Eq. ()2.1|) . except that the two root points 
must be labeled with i and j and the field point must be 
labeled with k. Consequently, yij\ k = fikfkj- Thus, Eq. 
(|4~T3l) yields 

Cij-kl = Vij-U + fVij-kl - - (likfy kjl £ + JitfVlj-k 

+JikfVki;i + JjifVu-k) ~ JikfkiJlj- (4.19) 

As for the function yij\ k t{jc), it is represented by the same 
diagrams as in Eq. \2.2\ , except that (a) the root points 
are labeled with i and j, while the field points are labeled 
with k and £, and (b) the first and second diagrams on 
the right-hand side of Eq. (|2.2[) . with their correspond- 
ing factors, actually become 2 and 4 labeled diagrams, 
respectively. In other words, each diagram is replaced by 
a symmetrized sum of topologically analogous diagrams, 
divided by its number, so that the symmetry property 
Vij\ki{ v ) — Vij\ek( r ) is preserved. We will refer to this 
process of generating the diagrams of the multicompo- 
nent case from those of the one-component case as sym- 
metrization. 
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Analogously to Eq. (|2.5I) . we can introduce the class at this point, it is at least clear from Fig. 12 of Ref. Hi] 



of approximations y^^ 2 \r). Therefore, Eqs. (14.1[) and 

,(Ai,A 2 ) 



(|47[) imply that B\^> is given by Eq. (|2?7)l . except 
that again each diagram is symmetrized to preserve the 
symmetry of under any permutation of indices. In 

the case of b\^^\ use of Eqs.|H8]and (f47T9|) shows that 
it is given by Eq. (I2.14[) . again with the symmetrization 
criterion. 

The properties (I2.16[) - (|2.18[) do not necessarily hold 
for individual labeled diagrams, but they do when again 
each diagram is replaced by its symmetrized sum. As a 
consequence, one has 



B 



(1,A 2 ) 



ij k£ , v 



3 _ B (*iAi) 



2 + Ai ijM > c 



A 2 = 



3Ai 
2 + Ai 



(4.20) 



This generalizes Eq. ([2~2"P]) . and hence Eq. pTT^|) . to the 
case of an arbitrary mixture. 



V. DISCUSSION 

Most of the integral equation theories are consistent 
with the exact correlation functions to first order in den- 
sity, Eq. (|2.1j) . so that they agree with the exact third 
virial coefficient. Therefore, the fourth virial coefficient 
is the earliest one that, not only differs from theory to 
theory, but even among different thermodynamic routes 
within the same theory. Thus, if one considers the two 
most studied liquid state theories (HNC and PY) and 
the three most standard thermodynamic routes (virial, 
compressibility, and energy) , there are in principle six al- 
ternative approximations for the fourth virial coefficient 
of a given interaction model. Actually, this number is 
reduced from six to five because the HNC theory belongs 
to the "exclusive" class of approximations that are ther- 
modynamically consistent with respect to the energy and 
virial routes,— ~— without being forced to do so. 

The main aim of this paper has been to prove that the 
number of independent predictions given by the HNC and 
PY approximations further reduces from five to four be- 
cause the fourth virial coefficient obtained from the virial 
route in the HNC theory is exactly three halves the value 
obtained from the compressibility route in the PY the- 
ory, for any interaction potential, any number of com- 
ponents, and any dimensionality. This result has been 
derived as a special case of a more general mathematical 
property described by Eqs. (|2.20p and (|2.21[) and graph- 
ically sketched by Fig. [TJ 

It is interesting to remark that, as a simple corollary of 
Eq. (|1.19l) , the Boyle-like temperature Tb 4 at which B4 = 
is the same in the PY-c and HNC-v approximations. 
Since both predictions cross at T = Ts 4 , none of them is 
closer to the exact value than the other one for the whole 
temperature range, as illustrated by Figs. HHU 

One may reasonably wonder if a relation similar to Eq. 
(|1.19j) extends to higher-order virial coefficients. While 
we are not in conditions of ascertaining this possibility 



that the ratio B™ c / Bf^ is not a constant. 

In conclusion, we expect that the simple relationship 
(|1.19j) can help to gain new insight into the general ther- 
modynamic inconsistency problem, as well as into the 
connection between the HNC and PY theories. From 
that point of view, the result (|1.19|) has a clear pedagog- 
ical value, especially taking into account the scarcity of 
exact and general results in statistical mechanics. On the 
other hand, from a more practical viewpoint, Eq. (|1.19p 
can be useful to test the correctness of analytical eval- 
uations and/or the accuracy of numerical computations 
of B4 in the HNC and PY frameworks, especially in the 
case of mixtures. 
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Appendix: Proof of Eqs. ([2716)) - ([2718)) 



In this Appendix we derive the identities (|2.16|) - (|2.18p . 
For simplicity, we will consider here zero-root diagrams, 
which are just V times the corresponding one-root dia- 
grams [cf. Eq. (|2.9j> ]. Also, we will introduce diagrams 
such that an arrow on a solid bond from point i to point 
j represents a term fijifij while an arrow on a dashed 
bond represents a term Vy/y. For instance, 




dl / d2 / d3 / d4/ 13 / 3 4/ 2 4/i4ri4-V 12 / 12 . 

(A.I) 



1. Proof of Eq. (|2~T6|) 

Let us start from 



dl / d2 / d3 / d4/ 13 / 34 / 2 4ri 2 • V12/12. 

(A.2) 

Now, taking into account that Vi 2 /i 2 = — V2/12, we can 
integrate by parts with the result 



(A.3) 
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where use has been made of the properties V2 • r±2 = 
—d and V2/24 = ^2if2A- Moreover, upon writing the 
second term on the right-hand side of Eq. (|A.3|) , we have 
exploited the invariance of the zero-root diagrams under 
rotation and reflection. Next, making Y2A = — 1*12 + ri3 + 
r-34, one has 



Integration by parts of the last term on the right-hand 
side yields 



From Eqs. (IA.4I) and (IA.5|) one gets 



1 1 


1 


1 1 




~3 





(A.5) 



(A.6) 



Substitution of Eq. (fA~6|) into Eq. (jA~3) gives Eq. (|27T6l) . 



2. Proof of Eq. (|2~TT)| 

Steps similar to those of Eq. (|A.3|) lead to 



(A.7) 



Next, setting ri 4 = ri 2 + r 2 4, one has 



(A. 4) Combining Eqs. (|A.7[l -(jA.9 |) we get 




(A.8) 



(A.9) 



(A.10) 



Finally, Eqs. (|A~7| and (|A~10t gives Eq. (|2Tf)) . 



3. Proof of Eq. (1238]) 

Again, integration by parts yields 




= -d 





Also, the equality r24 = — ri2 + 1*14 gives 






(A.ll) 



(A.12) 



Equation (I2.18[) directly follows from Eqs. (|A.11I) and 
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